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Supersymmetric Higgs Bosons in Weak Boson Fusion 
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We compute the complete supersymmetric next-to-leading order corrections to the production 
of a light Higgs boson in weak boson fusion. The size of the electroweak corrections is of similar 
order as the next-to-leading order corrections in the Standard Model. The supersymmetric QCD 
corrections turn out to be significantly smaller than their electroweak counterparts. These higher- 
order corrections are an important ingredient to a precision analysis of the (supersymmetric) Higgs 
sector at the LHC, either as a known correction factor or as a contribution to the theory error. 
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The main task of the LHC era is to understand elec- 
troweak symmetry breaking and the ultraviolet comple- 
tion of the Standard Model. According to electroweak 
precision data we expect to see a light Higgs boson, which 
should be embedded into a UV completion solving the 
hierarchy problem. A minimal realization of TeV-scale 
supersymmetry (MSSM) is a leading candidate for that. 
The most promising discovery channel for a light super- 
symmetric Higgs boson is the production in weak-boson 
fusion with a subsequent decay into tau leptons [l], 
There, a light supersymmetric Higgs scalar is guaranteed 
to appear over the entire MSSM parameter space 0. 

In the Standard Model as well as in the MSSM, many 
years of LHC running will be devoted to understand- 
ing the Higgs sector in detail, for example extracting the 
gauge and Yukawa couplings [J] . To meaningfully distin- 
guish between, for example, the Standard-Model and the 
MSSM Higgs sectors [5( we need to control the theory 
error on the LHC rates including higher-order effects. 
In the Standard Model the next-to-leading order QCD 
and electroweak corrections to weak-boson-fusion Higgs 
production are known to be fairly small [(J. In particu- 
lar, the QCD corrections are suppressed due to the color 
structure of the production process and the forward-jet 
topology. Also interference effects between Higgs produc- 
tion in weak-boson fusion and in gluon fusion with two 
additional jets are strongly suppressed Q. 

For either a comparison between the Standard-Model 
and the MSSM Higgs sectors or for a precision analysis 
of the MSSM Higgs sector these higher-order corrections 
have to be augmented by supersymmetric particle loops. 
In parallel to the Higgs searches, the LHC experiments 
will also search for direct signatures of new physics. If 
we should find such new states we can then predict their 
effects on the Higgs sector. If, for example, squarks and 
gluinos should be too heavy or the spectrum should be 
not favorable to precision MSSM analyses [g], we need to 
include their effects in the theory errors. Both cases re- 
quire a comprehensive calculation of the supersymmetric 
contributions to the weak-boson-fusion and gluon-fusion 
production processes 0]. 



Supersymmetry vs Standard Model — Compared to 
its Standard-Model counter part the leading-order pro- 
duction rate of a light supersymmetric Higgs scalar h 
includes an additional coupling factor sin(/3 — a). It is 
expressed in terms of the ratio of the vacuum expecta- 
tion values tan (J and the scalar mixing angle a from the 
supersymmetric two-Higgs-doublet model. For a given 
Higgs mass we can relate the tree-level MSSM production 
rate to the Standard-Model result via a simple re-scaling 
by sin 2 (/3 — a). For large pseudoscalar Higgs masses 
rriA ^ 200 GeV this factor is very close to unity. 

Including higher orders, there are additional contribu- 
tions from the supersymmetric particle spectrum: first, 
we take into account loops of supersymmetric partners. 
If we assume R parity, one-loop diagrams cannot mix su- 
persymmetric and Standard-Model particles, allowing for 
a diagram-by-diagram separation of the MSSM contri- 
butions. Secondly, the additional supersymmetric Higgs 
bosons with their Standard-Model type R charge appear 
in loops. Because the Standard-Model Higgs boson does 
not simply correspond to one supersymmetric Higgs bo- 
son we cannot separate the Standard-Model Feynman 
diagrams from the MSSM set. Instead, we first com- 
pute the MSSM-Higgs corrections and then subtract the 
Standard-Model Higgs loops, scaled by the tree- level cor- 
rection factor. 

Because of the large number of Feynman diagrams we 
compute the cross sections using the automated tool Had- 
Calc [10]. The Feynman diagrams and the amplitudes we 
generate with FevnArts/FormCalc II ll . The loop inte- 
grals we evaluate using LoopTools [12j . We assume min- 
imal flavor violation, because after taking into account 
all experimental and theoretical constraints, the effect of 
non-minimal flavor violation on LHC rates is small [l3| . 
We also assume a CP-conserving MSSM. 

Higgs-sector corrections — The mass of the light su- 
persymmetric Higgs boson mh is not a free parameter. 
At tree level it can be computed from and tan (3 
and is always smaller than mz- Higher-order correc- 
tions, dominated by the top Yukawa coupling term hf 
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effective theory 


Feynman 


diagrams 




Aa/a(ud — > udh) 


(0"a eft — CTfull)/cr 




Qeff full 


"off 


full 




effective theory 


\hhh 

*HHh 


r\ nrio f\ i no 

0.208 0.198 
—0.285 —0.275 


0.210 
-0.284 


0.210 
-0.279 


Hcff 

full 


-0.389 % 
-0.266 % 


-0.122 % 


Ajf hh 


-0.216 -0.219 


-0.220 


-0.257 




Feynman diagrams 


^hhh 


0.952 1.503 


0.950 


1.276 


Hcff 


-0.393 % 


-0.076 % 


«cff 


-0.1132 


-0.1158 


full 


-0.317 % 


m h 


109.8 GeV 


111.0 GeV 


Feynman diagrams, loop-improved Zfh 


m H 


391.5 GeV 


391.6 GeV 


OteS 

full 


-0.343 % 
-0.228 % 


-0.115 % 



Table I: Higgs self couplings for the parameter point SPSla 
following Ref. [JJ] (left) and Ref. [3 (right). The common 
factor — 3ernw / (2cw sw) is not included. 
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ush nih to values beyond the LEP2 
Phcnomenologically relevant 
studies therefore need to include quantum corrections to 
the Higgs mass and the Higgs potential. 

The challenge in including these higher-order correc- 
tions in our calculation is that we cannot simply shift the 
final-state Higgs mass. Already in the Standard Model 
the physical Higgs mass rtih is linked to the running quar- 
tic coupling X(Q) and top Yukawa h t {Q) flij ] 
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with Z(x) = 2^/xtan~ 1 (x~ 1/2 ) for x > 1/4. When in- 
cluding self-energy corrections to the Higgs mass, we 
should also correct the Higgs self couplings at the same 
order in perturbation theory. While this has to be 
done explicitely [20l | in a Feynman-diagrammatic ap- 
proach [3, EH it is automatically taken care of if we com- 
pute the quantum effects in the scalar potential [Til . 17 1. 

On the other hand, quantum corrections to the Higgs 
potential are usually computed for vanishing external 



momentum instead of p 



In the Feynman- 



diagrammatic approach it is straight forward to compute 
the Higgs self-energy diagrams with a finite momentum 
flow, while this is a challenge for the effective-potential 
method. 

Last but not least, when computing the effective scalar 
potential using renormalization-group techniques, it be- 
comes increasingly tedious to separate scales, like the 
heavy Higgs mass tua, the light stop mass and the 
gluino mass m~ g . 

Because the light supersymmetric Higgs boson is the 
lightest particle in the supersymmetric loops we expect 
the numerical effects of its mass and of the associated 
self coupling Xhhh to be non- negligible. In Tab. U we com- 
pare the values for the scalar self couplings in the effective 
theory approach (implemented from Ref. [T3| ) with those 
from the Feynman-diagrammatic FeynHiggs package [IH . 
Both approaches allow for an approximate computation 
introducing an effective mixing angle a e s from the scalar 



Table II: Schemes for computing the supersymmetric correc- 
tions to the VVh vertices at the hadronic level, for the leading 
partonic subprocess. 



Higgs mass matrix. In this case the trilinear Higgs cou- 
pling Xhhh is given by the Standard-Model coupling times 
the MSSM correction factor cos 2a e s sin(/3 + a e g). The 
results in this approximation should be equivalent in both 
schemes, because finite values of p 2 as well as corrections 
to the self couplings are skipped. In Tab. U we see that 
the a e ff approximations indeed agree very well. Between 
the full results, where FeynHiggs includes the fixed-order 
one-loop 0(hf) corrections to Xhhh [20| . we see the ex- 
pected small deviations. 

For LHC cross sections this comparison is compli- 
cated by the different final-state Higgs masses in the two 
schemes. The effect of all supersymmetric corrections 
to the VVh vertices (V = W, Z) and the Higgs wave- 
function renormalization for the dominant partonic sub- 
process ud — > udh we show in Tab. [Til The contributions 
from the Higgs sector and from supersymmetric particles 
cannot be separated because of their combined renormal- 
ization. The first four lines use a wave-function renormal- 
ization at one loop. In the a e ff approximation the cor- 
rections agree well, despite the different external Higgs 
masses. For the bottom lines we include the higher-order 
improved wave- function renormalization factors provided 
by FeynHiggs. 

For the numerical analysis in this letter we use the 
Feynman-diagrammatic approach with one-loop Z fac- 
tors, to allow for a proper renormalization-scale behav- 
ior. The difference of ~ 0.05 % between the different 
rate predictions is a lower limit on the remaining theory 
error from the MSSM Higgs sector. Note, however, that 
this theory error only applies if we strictly assume the 
minimal renormalizable supersymmetric Higgs sector. 

Supersymmetric particle corrections — Example one- 
loop diagrams appearing in the qq — ► qqh processes we 
show in Fig. Q] Our numerical results are based around 
the parameter point SPSla 2l|) for which there are two 
aspects we need to remember: supersymmetric particles 
with only electroweak charges have typical masses of 100- 
200 GeV, while squarks and gluinos range around 500- 
600 GeV, and tan/3 = 10 avoids large non-decoupling 
effects from down-type fermions. 



Figure 1: Left to right: Feynman diagrams contributing to strong vertex corrections, strong boxes, strong pentagons. 



As we can see in Tab. Mil all QCD corrections Act oc 
a 3 a s turn out to be surprisingly small. In the Stan- 
dard Model we know that from a QCD perspective we 
are essentially looking at two-sided non-interfering deep 
inelastic scattering. However, there are several mecha- 
nisms responsible for an even larger suppression in the 
supersymmctric case. 

Two tree-level vertices receive one-loop corrections, 
qqV and VVh, but only the first is corrected by 
squark/gluino loops. The qq'W coupling connects left- 
handed sfermions. Since the mixing between left and 
right-handed light-flavor squarks is proportional to the 
negligible quark Yukawa coupling, both external quarks 
are then left-handed, just as at tree level. This means 
that in the one-loop diagram (when closed with the 
Born diagram) the left-handed fermion trace cannot be 
connected through a gluino-mass insertion, because this 
would require a chirality flip. Instead of m^, the typical 
momentum scale in the numerator is ~ an order 

of magnitude below the gluino mass in the denominator. 

In the electroweak case, also the (typically lighter) 
charginos and neutralinos in the loop couple to the vector 
boson. This means that we can add a double mass inser- 
tion into the fermion line which can partly compensate 
for the heavy masses in the loop denominator. This effect 
leads to a relative enhancement of the electroweak over 
the QCD qqV vertex correction we observe in Tab. Mil 

For strongly interacting boxes, the qq'W and qqg cou- 
plings are the same for both diagrams shown in Fig. [TJ 
but the qqh coupling is proportional to T3 — Qs^,, i.e. 
around —1/3 for down squarks and +5/16 for up squarks. 
This leads to a cancellation by one order of magnitude, 



diagram Aa/a [%] 


diagram Aa/a [%] 


Aa ~ O(a) 


Aa ~ 0(a s ) 


self energies 0.199 




qqW + qqZ -0.392 
qqh -0.0260 
WWh + ZZh -0.329 


qqW + qqZ -0.0148 
qqh 0.00545 


box 0.0785 
pentagon 0.000522 


box -0.00518 
pentagon -0.000308 



sum of all Aa/a = -0.484 % 



Table III: Complete supersymmetric corrections to the pro- 
cess pp — > qqh by diagrams. Our parameter point SPSla has 
a tree-level rate of 706 fb. 



which could only be broken by different squark masses. 
Left-handed squarks, however, form a SU(2) doublet and 
are governed by the same soft-breaking terms, and the 
left -right mixing is negligible for light-flavor squarks. 
This argument does not hold for the sub-leading ZZ fu- 
sion, where we indeed find that the corrections turn out 
to be at a more natural level. 

In the Standard Model the color factor of a gluon ex- 
change between the two incoming quarks is proportional 
to the trace of the SU(3) generators and hence zero. 
The same is true for a pentagon gluino exchange between 
the incoming quarks, where the color trace is evaluated 
along quark/squark lines. In Fig.[TJ we show another su- 
persymmetric pentagon diagram with a squark exchange 
between the two incoming quarks. The VV- fusion is re- 
placed by a squark coupling to the Higgs, which gets rid 
of the color suppression. Such diagrams contribute for- 
mally at order 0(a^a 2 ), which is as large as the Born 
term (D(a 3 ). However, their kinematic properties are 
completely different from the vector-boson-fusion topol- 
ogy and the large loop masses further reduce their con- 
tribution to an altogether negligible level. 

Following all the above arguments the supersymmetric 
QCD corrections to weak-boson-fusion Higgs production 
are suppressed by a whole list of mechanisms, which ex- 
plain their at first sight surprising suppression even with 
respect to the electroweak corrections in Tab. Mil 

Looking beyond SPSla, we show the next-to-leading 
order corrections for the complete set of SPS parame- 
ter points [2l| in Tab. IIVI From the discussion above 
we do not expect the picture of electroweak vs. strong 
corrections to change significantly for any of them. Heav- 
ier supersymmetric spectra and different values of tan fi 
and of the trilinear couplings just scale the over-all size 
of the supersymmetric corrections. The relatively large 
corrections for the SPS5 parameter point are driven by a 
light top squark, while the largely decoupled spectrum in 
SPS9 leads to negligible MSSM effects. The typical size 
of the complete MSSM corrections is less or around 1 %. 

To study the behavior of the one-loop corrections with 
varying supersymmetric masses we start from the param- 
eter point SPSlb and run the universal gaugino mass 
toi/2 from 100 to 1000 GeV. In Fig. [2] we show the re- 
sult for a mT/2-dependent Higgsino mass parameter as 
well as for the fixed SPSlb value fx = 499 GeV. The 
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Aa/a [%} 





WWh + ZZh 


0(a) 


0(a s ) 


all 


SPSla 


-0.329 


-0.469 


-0.015 


-0.484 


SPSlb 


-0.162 


-0.229 


-0.006 


-0.235 


SPS2 


-0.147 


0.129 


-0.002 


-0.131 


SPS3 


-0.146 


-0.216 


-0.006 


-0.222 


SPS4 


-0.258 


-0.355 


-0.008 


-0.363 


SPS5 


-0.606 


-0.912 


-0.010 


-0.922 


SPS6 


-0.226 


-0.309 


-0.010 


-0.319 


SPS7 


-0.206 


-0.317 


-0.006 


-0.323 


SPS8 


-0.157 


-0.206 


-0.004 


-0.210 


SPS9 


-0.094 


-0.071 


-0.003 


-0.074 



Table IV: Complete MSSM corrections for all SPS parameter 
points [2lj . The vertex correction in the first column corre- 
sponds to TabHTl but including all partonic channels. 

corrections sharply drop with increasing m^, as we ap- 
proach the decoupling limit. Fixing \i means larger cor- 
rections for a light SUSY spectrum and a sharper drop 
for heavy masses. The maximum size for the corrections 
consistent with the LEP2 chargino limit we read off to 
be — 2 %. If we tune all weak-scale MSSM parameters to 
barely respect all LEP2 and Tevatron limits we find that 
the size of the supersymmetric corrections is bounded 
by —4%. Explicit non-decoupling effects in the bottom 
Yukawa only appear in this process at the two-loop level, 
which means all curves in Fig. [3 decouple smoothly for 
increasing masses. Consistent with our previous discus- 
sion the 0{a s ) corrections are negligible over the entire 
parameter range. 

Outlook — In the light of a possible precision analysis 
of the Standard-Model and MSSM Higgs sector at the 
LHC we have analyzed the size of the supersymmetric 
one-loop corrections to the weak-boson-fusion production 
process qq — > qqh. 

The appearance of all supersymmetric neutral Higgs 
bosons in the loops required us to study the impact of 
different methods of describing higher-order effects on 
Higgs masses and the Higgs potential. We find that the 
corrections from the Higgs sector are at the per-cent level, 
with a remaining uncertainty of below 0.1 % due to these 
calculational approaches — simply reflecting unknown 
higher-order corrections. 

The supersymmetric one-loop QCD corrections are not 
only suppressed to a typical NNLO level, but turn out to 
be negligible. This is due to a variety of effects, based on 
the color structure, the supersymmetric coupling struc- 
ture, or the kinematics of the process. The complete set 
of electroweak loop diagrams contributes at the per-cent 
level, as is expected for massive 0(a) corrections. 

In total, the supersymmetric one-loop corrections to 
Higgs production via vector-boson fusion can be up to 
4 % for parameter points allowed by direct SUSY searches 
and are typically at or below 1 %. Their sign is in general 
negative. This result should serve as a solid basis for a 
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Figure 2: Relative next-to-leading order corrections as a func- 
tion of m-t/2 for varying and for fixed fi. For the latter we 
show the strong corrections independently. The vertical lines 
indicate the chargino mass limit from LEP2 and the reference 
point SPSlb. 



precision analysis of the supersymmetric Higgs sector at 
the LHC. 
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